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ABSTRACT 

Using samples drawn from the Sloan Digital Sky Survey, we study the relationship between local 
galaxy density and the properties of galaxies on the red sequence. After removing the mean dependence 
of average overdensity (or "environment" ) on color and luminosity, wc find that there remains a strong 
residual trend between luminosity-weighted mean stellar age and environment, such that galaxies 
with older stellar populations favor regions of higher overdensity relative to galaxies of like color and 
luminosity (and hence of like stellar mass). Even when excluding galaxies with recent star-formation 
activity (i.e., younger mean stellar ages) from the sample, we still find a highly significant correlation 
between stellar age and environment at fixed stellar mass. This residual age-density relation provides 
direct evidence for an assembly bias on the red sequence such that galaxies in higher-density regions 
formed earlier than galaxies of similar mass in lower-density environments. Wc discuss these results 
in the context of the age-metallicity degeneracy and in comparison to previous studies at low and 
intermediate rcdshift. Finally, we consider the potential role of assembly bias in explaining recent 
results regarding the evolution of post-starburst (or post-quenching) galaxies and the environmental 
dependence of the type la supernova rate. 

Subject headings: galaxics:statistics, galaxies: abundances, galaxies:stellar content, large-scale struc- 
ture of universe 



1. INTRODUCTION 

Over the past decade or more, the formation his- 
tories and clustering properties of dark matter halos, 
within a Lambda cold dark matter (ACDM) cosmol- 
ogy, have been studied in detail using both numerical 
simulations and analytical analyses. These theoretical 
studies have repeatedly shown t hat more massive halos 
are more strongly c l ustered (e.g., lEfstathiou et al.lll988l; 
Cole fc Kaiserfl989t[Mo~fc Whitdll996l : lSheth fc Tormenl 



19991 : iWetzel et all 120071 ). a result that is supported 



by the observed clustering of nearby galaxy groups 
dPadilla et al.|[200l lYang et al.ll200l iBerlind et afll2006l : 
iWang et al.ll2008T L Furthermore, recent dark matter sim- 
ulations indicate that the clustering of halos of a given 
mass depends on their assembly history, where halos 
that assembled their mass earlier in the history of the 
Universe are more strongly clus tered than halos of the 
same mass that formed later (jSheth fc Tormenl 12004 : 
iGao et aT]|2005f): this effect is com monly referred to as 
assembly bias (jCroton eFaD[2003). 

The concept of assembly bias could be applicable to the 
formation of galaxies in addition to dark matter halos, 
such that galaxies that assembled their stellar mass ear- 
lier would be more strongly clustered today than galaxies 
of like mass and younger stellar populations. Such a bias 
in the galaxy population has typically been explored in 
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two distinct manners: [1] by studying galaxy samples 
spanning a range of rcdshifts and directly observing the 
evolution in the population as a function of mass and 
clustering properties, and [2] by applying a more archae- 
ological or paleontological approach, where the stellar 
populations of nearby galaxies are studied in detail, with 
the goal of disentangling the typical evolutionary history 
in different environments. 

Various observational studies have investigated the for- 
mation of early-type or red-sequence galaxies from both 
perspectives. For instance, analyses of sa mples drawn 
from the DEEP2 Galaxy Red shift Survey dDavis et all 
120031 iNewman et all 12009ft by iBu ndv et all (|2006|) and 
bv iCooper et al l (|2006l . I2007T) showed that the popu- 
lation of galaxies on the red sequence at z < 1 was 
preferentially built-up in overdense environments, so 
that at a given mass red galaxies in dense environs 
have typically ceased their star formation (i.e., assem- 
bled their mass) earlier than those in less dense regions. 
This picture is supported by analyses of the fundamen- 
tal p lane (FP, iDiorgovski fc Davisl 119871 : iDressler et~aTl 
119871 ) at similar redshifts, where it has been found that 
early-type galaxies in high-density regions have reached 
the FP more quickly than those in lower-density en- 
virons (Ivan Dokkum et all 120011: iGebhardt et all (20031 : 
iTreu et al.ll2005HMoran et al.ll2005ft . ~ 

In contrast to the relatively coherent picture derived 
from studies of early-type or red-sequence galaxies at 
intermediate redshift, detailed studies of the stellar pop- 
ulations in nearby ellipticals have yielded mixed results. 
A recent comparison of the stellar ages of early-type sys- 
tems in t he Coma cluster ve rsus a corresponding field 
sample by iTrager et al.l (|2008f ) found no difference in as- 
sembly history between the two environment regimes. 
On the other hand, a similar analysis of elliptical and 
lenticular galaxies in the Virgo and Coma clusters by 
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iThomas et al.l (|2005() found a distinct relationship be- 
tween mean stellar age and the environment in which 
a galaxy resides, such that cluster early-types are older 
than their field counterparts. 

In this paper, we take a different approach towards 
studying the topic of assembly bias, with hopes of rec- 
onciling the results from studies at intermcdiatc-rcdshift 
with those focusing on local samples. Many previous ef- 
forts at low redshift relied on comparing cluster and field 
populations and/or focused on high quality observations 
of a relatively small number of galaxies and/or partially 
controlled for morphology while ignoring inter-sample 
variations in properties such as rest-frame color. In con- 
trast, we utilize somewhat lower-quality data for a much 
larger sample of g alaxies drawn fro m the Sloan Digital 
Sky Survey (SDSS. lYork et al.ll2000h to study the assem- 
bly histories of thousands of red galaxies in a broad and 
continuous range of environments, controlling for correla- 
tions between environment and color, morphology, mass, 
etc. In the following Section ($2]), we describe the galaxy 
sample under study, including measurements of galaxy 
environments and stellar ages. Our results regarding the 
relationship between stellar mass, age, and environment 
are then presented in Ej3j followed by further analysis, 
discussion, and a summary in SJSJ and respec- 
tively. Throughout, we employ a ACMD cosmology with 
w = —1, Q m = 0.3, Q\ = 0.7, and a Hubble parameter 
of Hq = 100 h km s _1 Mpc -1 , unless otherwise noted. 

2. DATA 

With spectra and multi-band photometry for more 
than 500, 000 galaxies, the SPS S Data Release 4 (DR4, 
lAdelman-McCarthv et aLll2006f ) enables the local density 
of galaxies at z < 0.2 to be measured over nearly one 
quarter of the sky. We select a parent sample of 372, 688 
galaxies from the SDSS DR4, as drawn from the New 
York University Value-Add ed Galaxy Catalog (NYU- 
VAGC lBlanton et al.l2005bf ). This sample is drawn from 
the NYU-VAGC Large-Scale Structure (LSS) catalog, 
which provides a relatively uniform spatial selection func- 
tion by combining the SDSS imaging, target, and tiling 
masks while tracking the flux limit and completeness as 
a function of angular position. Our parent sample is lim- 
ited to the redshift range 0.01 < z < 0.3 and to SDSS 
fiber plates for which the redshift success rate for targets 
in the main spectroscopic survey is 80 percent or greater. 

Rest-frame colors, absolute magnitudes, and stel- 
lar masses are computed using the KCORRECT K- 
correc tion code (version v4_l_4) of IBlanton fc Roweii 
(|2007l see also IBlanton et al.l l2003al) . The rest-frame 
quantities for the SDSS samples are derived from the 
SDSS apparent ugriz model magnitudes, where all mag- 
nitudes in this pap er are calibrated to the AB system 
(|Oke fc Gunnl 119831 ). Stellar masses are computed us- 
ing templat e spectral energy di s tribut ions (SEDs) based 
on those of IBruzual fc Charlotl (|2003l ). The best-fitting 
SED, given the observed ugriz photometry and spectro- 
scopic redshift, is used to directly co mpute the stellar 
mass-to-light ratio (M„»/L), assuming a lChabrierl (|2003f ) 
initial mass function. 

As shown by many previous studies at low and inter- 
mediate redshift (e.g., IStrateva et al.ll200lHBaldrv et al.1 
l200i [Mret~aT1l2004bl: IWillmer et al.N2006D . and as il- 
lustrated in Figure [TJ the distribution of galaxies in 



color-magnitude space is bimodal, with a relatively tight 
red sequence and a more diffuse blue cloud. To iso- 
late the red-sequence population, we use the following 
magnitude-dependent cut: 

g-r = -0.02667- M r + 0.11333. (1) 

This division in g — r color is shown in Figure [1] as the 
dashed red line. 
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Fig. 1. — We plot the rest-frame g — r versus M r color-magnitude 
distribution for SDSS galaxies in the spectroscopic sample within 
the redshift range 0.05 < z < 0.15. Due to the large number of 
galaxies in the sample, we plot contours (rather than individual 
points) corresponding to 25, 50, 200, 500, 1000, and 2000 galaxies 
per bin of A(g - r) = 0.05 and AM r = 0.1. The dashed red 
horizontal line shows the division between the red sequence and 
the blue cloud as given in Equation [T] 

To quantify stellar mctallicities (Z) and ages (£) 
al ong the red sequence , we employ the measurements 
of iGallazzi et al.l (|2005f ). which are based on model fits 
to spectral absorption features in the SDSS DR4 spec - 
tra. The model spectra utilized by IGallazzi et al.l (|2005f ) 
are derived from t he po pulation synthesis models of 
IBruzual fc Charlotl (|200l and span a broad range of 
star-formation histories. The models are simultaneously 
fit to a minimum set of metal- and age-sensitive spectral 
indices (D4000, H/3, US a + [Mg 2 Fe], and [MgFe], 

iFaber et al.|[l98 5: Worthcy 1994), yielding measurements 
of age and metallicity with typical uncertainties on the 
order of a t < 0.15 dex and oz < 0.3 dex, respectively. 

While both the age and metallicity measurements are 
sensitive to the signa l-to-noise ratio of th e galaxy spec- 
trum (see Table 1 of IGallazzi eTaLl [20051 ) . we select our 
sample down to a relatively-faint magnitude limit of 
r = 17.5. We investigate the sensitivity of our results 
to this apparent r-band magnitude limit in fj4j Among 
the 114,916 galaxies with r < 17.5 on the red sequence 
at 0.05 < z < 0.15, more than 93% have robust age 
and metallicity measures. The ~ 7% excluded from our 
sample populate specific (and narrow) redshift windows, 
where critical age- or mctallicity-constraining spectral 
features (e.g., H/3 or H<5) fall on sets of bright sky emis- 
sion lines. Thus, the galaxies with (or without) age and 
metallicity information are not biased with regard to as- 
sembly history (i.e., not biased in terms of stellar mass, 
morphology, etc.). 

In addition to the stellar mass estimates discussed 
ab ove, we also employ the stellar mass measurements 
of IGallazzi et al.l (|2005l ) , which differ from those com- 



Galaxy Assembly Bias on the Red Sequence 



3 



puted using KCORRECT in that they do not rely on 
fitting SEDs to t he SPSS photometry ; instead, the stel- 
lar masses from iGallazzi et al.l (|2005[ ) are a product of 
the same model fits to the SDSS spectra that are used 
to constrain stellar age and metallicity, with the SDSS 
model apparent magnitudes used only to normalize the 
resulting stellar mass estimates. Thus, these additional 
stellar mass values are complementary to those derived 
from KCORRECT while also being congruous with the 
age and metallicity values employed in ou r analy sis. Note 
that the s tellar masses of IGallazzi et al.1 ()2005l ) utilize a 
iChabrierl (|2003f ) initial mass function, the same as that 
used by KCORRECT. 

In the SDSS, we estimate the local galaxy overden- 
sity, or "environment" (as we will refer to it in this pa- 
per), using measurements of the projected fifth-nearest- 
neighbor surface density (£5) about each galaxy, where 
the surface density depends on the projected distance 
to the fifth-nearest neighbor, -Dp, 5, as £5 = 5/(7rD^ 5 ). 
Over quasi-linear regimes, the mass density and galaxy 
density shoul d simply differ by a factor of the galaxy bias 
(|Kaiserlll987f ). In computing £5, only objects within a 
velocity window of ±1500 km s _1 are counted, to ex- 
clude foreground and background galaxies along the line 
of sight. 

To correct for the redshift dependence of the sampling 
rate of the SDSS, each surface density is divided by the 
median £5 for galaxies within a window of Az = 0.02 
centered on the redshift of each galaxy; this converts the 
£5 values into measures of overdensity relative to the me- 
dian density (given by the notation I + S5 herein) and ef- 
fectively a ccounts for the reds hift variations in the selec- 
tion rate (jCooper et al. 200jl). We restrict our analyses 
to the redshift range 0.05 < z < 0.15, avoiding the low- 
and high-redshift tails of the SDSS galaxy redshift distri- 
bution where the variations in the survey selection rate 
are greatest. Finally, to minimize the effects of edges and 
holes in the SDSS survey geometry, we exclude all galax- 
ies within 1 h~ x Mpc (comoving) of a survey boundary, 
reducing our sample size to 102, 601 red galaxies within 
the redshift range 0.05 < z < 0.15. Note that all galaxies 
included in the NYU-VAGC LSS sample (independent of 
color and down to the r < 17.77 spectroscopic limit) are 
used to trace the local galaxy environment. 

3. RESULTS 

A wide array of galaxy properties at low and intermedi- 
ate redshift have been shown to correlate with environ- 
ment. For example, at z < 1, there is a strong color- 
density relation such that red galaxies tend to reside 
in hi gher density environs than their blue counterparts 
(e.g..lHogg eTaLl[200l iBalogh et al1[200l iCooper et all 
120061 ). Furthermore, within the red galaxy population, 
there is also a significant luminosity dependence to the 
clustering of galaxies, with the most luminous members 
of the red sequen ce residing in the most massive dark 
matter halos (e g iNorberg et al ll2002l : |Hogg et al 1120041 : 
iZehavi et aDl2005l : ICooper et al.lbo'oa ICoil et al.ll2006P ). 
This relationship between environment, color, and lumi- 
nosity is such that at fixed color and luminosity no sig- 
nificant trend is observed between local galaxy density 
and a variety of other galaxy properties s uch as surface 
brightness, Sersic index, or stellar mass (Blanton et al. 
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Fig. 2. — (a) For the population of galaxies on the red sequence, 
we show the mean galaxy overdensity, < log 10 (l + <%) >, as a 
function of rest-frame color, g — r, and absolute magnitude, M r , 
computed in a sliding box of height A(g — r) = 0.03 and width 
AM r = 0.11. The size and shape of the box are illustrated in the 
upper left-hand corner of the plot, (b) We plot the mean resid- 
ual environment, < A5 >. as a function of color and magnitude, 
computed in the same sliding box. 

However, a recent analysis of the relationship between 
environment and the luminosities, rest-frame colors, stel- 
lar masses, and gas-phase metallicities of nearby star- 
forming galaxies selected from the SDSS bv lCooper et al.l 
(|2008a| ) shows that after removing the mean dependence 
of environment on color and luminosity there remains a 
highly-significant correlation between environment and 
oxygen abundance. That is, they find that within the 
star-forming population (i.e., on the blue cloud) there 
is a relationship between local galaxy density and gas- 
phase metallicity separate from that observed between 
galaxy densit y and color and lumino sity (or stellar mass). 
Equivalently, ICooper et alJ (|2008al ) show that at fixed 
stellar mass and star-formation rate, star-forming galax- 
ies that reside in more overdense environments tend to 
have more metal-rich gas in their HII regions. 

In this paper, we apply the same techniques to study 
the relationship between mean stellar age and environ- 
ment a t fixed stellar mass a long the red sequence. Fol- 
lowing ICooper et al.l (|2008af) , we first measure and re- 
move (subtract) the mean dependence of environment 
on color and luminosity along the red sequence. Fig- 
ure [2^ shows the mean overdensity as a function of 
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rest-frame g — r color and r-band absolute magnitude, 
< log 10 (l + S 5 )[g — r,M r ] >, for our sample of red 
galaxies. In addition to a relatively weak color-density 
trend (when restricting to red galaxies), we find a strong 
luminosity-density relation, with the mean overdensity 
increasing at the under- and ov er-luminous ends o f the 
red sequence, as ob served by iHogg et all (|2003l) and 
iBlanton et all (|2005af ) . To remove this relationship of en- 
vironment to color and luminosity, we subtract the mean 
overdensity at the color and luminosity of each galaxy 
from the measured overdensity: 

A 5 =log 10 (l + 5 5 )- < log 10 (l + <S 5 ) [5 -r,M r ] >, (2) 

where the distribution of mean environment with color 
and absolute magnitude, < log 10 (l + <5s)[g — r, M r ] >, is 
median smoothed on A(<? — r) = 0.09 and AM r = 0.33 
scales prior to subtraction. Using rest-frame u — r color 
in lieu of g — r creates a greater dynamic range in color 
along the red sequence, but precipitates no change in our 
results regarding the relationships between environment 
and age, mctallicity, and stellar mass. 

The "residual" environment, A5, quantifies the over- 
density about a galaxy relative to galaxies of similar 
color and luminosity, where values of A5 greater than 
zero correspond to galaxies in environments more over- 
dense than that of the typical galaxy with like color and 
luminosity (i.e., with like star-formation history). Fig- 
ure O) shows the dependence of mean A5 on color and 
luminosity; no significant color or luminosity dependence 
is evident, with deviations from < A 5 > = being small. 
Fo r more details and di scussion of this technique, refer 
to iCooper et al.l (|2008af) . Finally, note that the "resid- 
ual" overdensity, A5 , is only a small perturbation to the 
"absolute" overdensity, log 10 (l + 65), such that the two 
environment statistics are still strongly correlated with 
each other for individual objects. 

In Figure [3J we show the average dependence of abso- 
lute environment, log 10 (l+#5), on luminosity, rest-frame 
color, stellar mass, stellar age, and stellar metallicity for 
galaxies on the red sequence. As seen in many previous 
studies, wc find a strong relationship between luminos- 
ity and local galaxy density, driven by luminous galaxies 
at the centers of groups and clusters and fainter satel- 
lite galaxies sitting in the outskirts of overdense regions 
(jBerlind et al.l[200a l. The color-density relation is rel- 
atively weak over the narrow range of rest-frame colors 
spanned by the red sequence. The small measured varia- 
tions in the average overdensity with color are likely asso- 
ciated with the projection of the luminosity-density rela- 
tion that results from the intrinsic tilt of the red sequence 
towards redder colors at brighter luminosities. In addi- 
tion, the observed fluctuations in the color-density rela- 
tion may be caused by interloping disk galaxies, which 
appear on the r ed sequence due to dust and inclina- 
tion effects (e.g . iWeiner et all 120051 iDriver et all 120061 ; 
lLotz et al.|[2008l ; iMaller et al.ll2009fk away from the ridge 
of the red sequence, these interlopers have a stronger im- 
pact on the average environment, due to their greater rel- 
ative numbers. Finally, the mass-environment relation- 
ship within our sample closely follows that between lu- 
minosity and environment, as expected given the strong 
correlation betwe en the combination of g — r and M r with 
stellar mass fe.g.. lBell &: de Jom3l200lHKauffmann et al.1 



[200llCooper et al.ll2008bh . 

In addition to strong luminosity- and mass- 
environment trends along the red sequence, we find a sig- 
nificant relationship between mean absolute environment 
and luminosity- weighted mean stellar age, where galax- 
ies with older stellar populations favor overdense regions 
(cf. Figure [3]). With a typical uncertainty of <r t < 0.15 
dex for each indiyidual age measurement (see Table 1 
of iGallazzi et alj [20051 ) . the residual age-density rela- 
tion spans multiple independent bins, while also covering 
nearly 0.4 dex in < A5 >. 6 

We also find a significant trend between stellar metal- 
licity, Z, and absolute environment along the red se- 
quence, such that more metal-rich galaxies favor higher- 
density regions. However, it should be noted that there 
is a non-negligible correlation between the binned data 
points in Fig. (3[ given the relatively large uncertainty in 
the typical metallicity measurement, az ~ 0.1-0.5 dex. 
Furthermore, the vast majority (> 75%) of the red galaxy 
population has measured mctallicitics of log lo (Z/Z0) > 
—0.1, so the apparent strength of the metallicity-density 
relation is driven in part by a relatively small num- 
ber of galaxies in the sample and by those systems 
with the greatest uncertainty in Z — the <tz values are 
typically larger at low metallicities, with a significant 
tail to a z ~ 0.5 dex. Altogether, the strength of the 
metallicity-density relation is considerably weaker than 
that of the luminosity-, stellar mass-, or age-density re- 
lations, though still highly statistically significant. 

By studying the dependence of residual environment, 
A 5 , on each galaxy characteristic (e.g., M*, t, and Z), 
we can determine whether any particular property ex- 
hibits an excess correspondence with local galaxy den- 
sity beyond that explained by the color-luminosity- 
environment relation. In Figure 21 we show the average 
dependence of mean residual environment on each galaxy 
property shown in Figure [3J (color, luminosity, stellar 
mass, age, and metallicity). By construction, there is no 
significant correlation between average A5 and luminos- 
ity or rest-frame color (cf. Table [T]). Similarly, we find no 
dependence of mean residual environment on stellar mass 
within the red-sequence population for both sets of stel- 
lar mass measurements described in The previously- 
noted tight relationship between stellar mass and rest- 
frame color and luminosity mass makes this result more 
or less expected. 

While wc find no significant correlation between resid- 
ual environment and color, luminosity, or stellar mass 
among red galaxies in the SDSS, we do see a strong re- 
lationship between A5 and stellar age, such that galax- 
ies with greater stellar ages tend to reside in regions of 
higher density relative to other galaxies of like color and 
luminosity (i.e., of like stellar mass). In other words, we 
find that at a given stellar mass red galaxies with older 
stellar populations favor overdense environments relative 
to their younger counterparts. 

Likewise, we find a significant correlation between stel- 
lar metallicity, Z, and residual environment; this resid- 
ual metallicity-density trend is such that more metal- 
rich galaxies favor more overdense regions than galaxies 

6 The median uncertainties in the age and metallicity measure- 
ments for the r < 17.5 galaxy sample are at = 0.11 dex and 
az = 0.13 dex, respectively. 
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Fig. 3. — The dependence of mean overdensity, < log 10 (l + 65) >, (i.e., "absolute" environment) on the absolute magnitude, rest-frame 
color, stellar mass, stellar age, and stellar metallicity of red-sequence galaxies. The points and error bars give the mean environments and 
la uncertainties in the means computed in distinct bins. In the top right panel, we show the results based on two independent stellar mass 
estimates, th ose determined solely from fits to the broad-band photometry using KCORRECT and those based on model fits to the SDSS 
spectra from Gallazzi ct al. (2005). As seen in previous studies, we find strong luminosity-density and mass-density relations along the 
red sequence. We also observe relatively strong correlations between average environment and other properties of the stellar population, 
such as age and metallicity. 



of like color and luminosity. The residual age-density 
and mctallicity-density relations are quite strong, span- 
ning roughly 0.2 dex and > 0.1 dex in A5, respectively, 
or more than half the range in overdensity spanned by 
the absolute age-density and metallicity-density rela- 
tions. However, as stated previously, the strength of the 
metallicity-density trend is driven by a relatively small 
number of galaxies at low metallicities (Log 10 (Z/Zo) < 
—0.1); in the more metal-rich regime, where the bulk of 
our sample resides, the residual metallicity-density rela- 
tion is relatively weak. 

4. POTENTIAL SYSTEMATIC EFFECTS 

In we presented the relationships between (both 
absolute and residual) environment and an assortment 
of galaxy properties, including 1 — band luminosity, rest- 
frame g — r color, stellar mass, luminosity- weighted mean 
stellar age, and stellar metallicity for a sample of red 
galaxies with r < 17.5. While we find no correlation be- 
tween mean residual environment and stellar mass, we 
do observe significant relationships between < A5 > and 
stellar age and metallicity. In the remainder of this sec- 
tion, we discuss possible systematic effects that could 
bias these results. 

4.1. Sample Selection and Aperture Effects 

In contrast to other studies of the relationship between 
environment and stellar age or metallicity for early-type 



TABLE 1 
Statistical Strengths of Relations 
with A5 



Galaxy Property 



slope 



Al r 
g — r 

log 10 (M./M e ) 

l°gio(*/yr) 
log lo (Z/Z ) 



^slopc 



-0.006 
-0.001 
0.009 
0.069 
0.039 



0.001 
0.046 
0.011 
0.333 
0.135 



0.002 
0.109 
0.005 
0.042 
0.017 



Note. - We quantify the correla- 
tions between various galaxy properties and 
residual density (A5) for the sample of red- 
sequence galaxies with r < 17.5. Here, we 
list the Spearman ranked correlation coef- 
ficient (p) as well as the measured slope 
(d < A5 > /dx) and its uncertainty derived 
from linear fits to the binned data points 
in Figure [4] The statistical strengths given 
in the table for the stellar mass-As relation 
correspond to the stellar mass values pro- 
duced by the KCORRECT code. 



galaxies in the local Universe (e.g., iTrager et al.1 120081 ). 
our work relies upon estimates of age and metallicity 
derived from relatively low-quality (i.e., low signal-to- 
noise ratio) spectra, which translates into measurements 
with greater uncertainties. It is only by constructing 
large samples that we are able to overcome these greater 
uncertainties and conduct detailed analyses of the mean 
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Fig. 4. — The dependence of mean residual environment, < A5 >, on absolute magnitude, rest-frame color, stellar mass, stellar age, 
and stellar mctallicity. The points and error bars correspond to the mean environments and la uncertainties in the means computed in 
distinct bins of each galaxy property. In the top right panel, we show the results based on two independent stellar mass estimates, those 
determined so l ely fro m fits to the broad-band photometry using KCORRECT and those based on model fits to the SDSS spectra from 
Gallazzi ct al. (2005). While wc find no significant residual trend between environment and color, luminosity, or stellar mass, we observe a 
strong dependence of residual environment on stellar age and mctallicity. 



relations with environment. 

Given the larger uncertainty in each individual age 
and metallicity measurement, it is important that we 
clearly understand any potential bias in our sample se- 
lection. The e r rors i n the estimates of both Z and t by 
IGallazzi et al.l (|2005l ) depend on the signal-to-noise ra- 
tio of the SDSS spectra, with metallicity measurements 
showing greater sen sitivity to the quality of the data. 
IGallazzi et al.l (|2005l ) find that the uncertainty in the es- 
timated age shows less variation with the mean signal- 
to-noise ratio per pixel of the SDSS spectrum than the 
uncertainty in t he metallicity constra int (cf. Section 2.4.1 
and Figure 6 of IGallazzi et al.|[20 05'). This is attributed 
to the fact that stellar age is largely constrained by 
the 4000A break and the Balmer absorption lines, while 
metallicity estimation relies heavily on Mg and/or Fe 
features, which are more difficult to measure at a given 
signal-to-noise ratio, especially in low-mass/metal-poor 
systems where these features are commonly weak. 

To test the dependence of our results on the uncer- 
tainty in the age and metallicity estimates (i.e., on the 
quality of the spectra), we vary the magnitude limit uti- 
lized in the sample selection, defining brighter subsam- 
ples for which the typical spectrum has a higher signal- 
to-noise ratio (cf. $2]). In addition to utilizing a r < 17.5 
limit, we study samples restricted to have r < 17 or 
r < 16.5. Across all of these samples, our general re- 
sults regarding the relationship between stellar age and 



environment remain unchanged; we still find a highly 
significant residual age density relation in each of the 
magnitude-limited samples (cf. Figure[5]), while no corre- 
lation is found between average A5 and color, luminosity, 
or stellar mass. Furthermore, the strength of the resid- 
ual age-density relation remains unchanged, even when 
limiting the sample to only those galaxies with higher- 
quality spectra (cf. Figure [5] and Table [2|) . 

While our results regarding the age-environment rela- 
tion are very robust to variations in the apparent mag- 
nitude limit of the sample, the relationship between stel- 
lar metallicity and environment becomes increasingly less 
significant as the sample is restricted to brighter (and 
therefore preferentially more luminous and more mas- 
sive) systems. In particular, the range in stellar mctal- 
licity probed by the r < 17 and r < 16.5 samples is 
reduced to \og 10 (Z/ZQ) > —0.75 and the uncertainties 
in the mean residual environment, < A5 >, increase such 
that the resulting correlation between average residual 
environment and stellar metallicity is only statistically 
significant within the relatively metal-rich regime (i.e., 
\o Sw (Z/Zq) > -0.35). 

In addition to concerns regarding the impact of low- 
quality data, our results could be biased by the in- 
clusion of dusty, edge-on star-forming galaxies in our 
analysis, which would skew the galaxy sample towards 
lower stellar ages in low density environments. Dusty 
star-forming systems may also be a concern in the 
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Fig. 5. — The dependence of mean residual environment, < A5 > 
on stellar mctallicity (top) and stellar age (bottom) , for a variety of 
sample selections. The black circles connected by the solid black 
line correspond to the main r < 17.5 red galaxy sample. The 
relations for three more-selective samples are ovcrplottcd: a r < 
16.5 magnitude-limited sample corresponds to the red diamonds 
connected by the dotted red line, a r < 16.5, n r > 2.5 sample 
corresponds to the blue squares connected by the dashed blue line, 
and a r < 17, n r > 2.5, b/a > 0.7 sample corresponds to the green 
triangles connected by the dashed-dotted green line. Even when 
including strong morphological, color, and brightness limits, the 
residual age— density relation persists. 

higher-density regime; recent results from the Space 
Telescope A901/902 Galaxy Evolution Survey (STAGES 
IGrav et al.ll200a IGallazzi et all 120091 : IWolf et all 120091) 
and the Local Clu ster Substructure Survey (LoCuSS, 
IHaines et al.l 12009! ) have suggested that obscured star 
formation may be prevalent in higher-density environ- 
ments such as the outskirts of clusters. As highlighted 
earlier in this paper, analyses of galaxy samples at 
intermediate and low redshift have shown that some 
galaxies on the red sequence are not truly "red and 
dead" , early-type systems, but are dust y (i.e., highly red- 
dened) star-formin g disk galaxies (e.g.. lBell et al.ll2004al : 
iWeiner et al.l |2005[ ) . This interloper population com- 
prises an increasingly grea ter fraction of the red sequence 
at low er luminosities (e.g.. lLotz~c t al. 2008: Maile r et al.l 
[20091) . so restricting our analysis to brighter magnitude 
limits, as discussed above, would be expected to reduce 
the impact of dusty disk galaxies. 

However, restricting our sample according to morphol- 
ogy provides a more robust means for minimizing the 
contribution of dusty disk galaxies to the measured age- 
density r^ki ; dOT^_Utili^mg_the r-band Sersic indices, 
n r , of iBlanton et al.l (|2003bf ). we select from our red 
galaxy samples those sources with n r > 2.5. For this 
color- and morphology-selected sample, the residual re- 
lationship between environment and stellar age remains 
strong. Even when restricting our analysis to r < 16.5 



and n r > 3, we still find no significant change in the 
strength of the residual age-density relation, though the 
uncertainties are greater due to the smaller sample size. 

To improve morphological discrimination further, we 
also employ measurements of the r-band isophotal ma- 
jor and minor axes from the SDSS imaging to remove 
highly-inclined disk galaxies. Restricting to sources with 
an axis ratio b/a > 0.7, with or without a joint morpho- 
logical selection (n r > 2.5), we continue to find a strong 
correlation between mean stellar age and residual envi- 
ronment (cf. Figure [5] and Table [2]). Furthermore, we 
also observe no significant dependence of mean A5 on n r 
or b/a, indicating that the observed correlation between 
stellar age and residual environment is not reflecting an 
unaccounted-for morphological dependence. 

While the age-density relation is quite robust to addi- 
tional selection criteria, the relationship between resid- 
ual environment and stellar metallicity becomes increas- 
ingly weaker and less significant as the galaxy sample is 
restricted according to color, apparent magnitude, and 
morphology. As shown in Figure [5l when limiting our 
sample to those galaxies with r < 17, n r > 2.5, and 
b/a > 0.7, we find weaker (though still statistically sig- 
nificant) evidence of a residual metallicity-density rela- 
tion, with an average residual environment of < A5 > ~ 
across the entire range of metallicities probed. 

A second way to test the robustness of our results 
to contamination on the red sequence is to make our 
magnitude-dependent color cut more restrictive, which 
will reduce the number of galaxies in the "green valley" 
or the red extremes of the blue cloud that are included 
in our sample. In addition to the brighter r-band mag- 
nitude limit of r < 16.5 and the Sersic selection cut of 
n r > 3, we shift our color selection, as given in Equation 
[TJ redward by g — r = 0.05. For this restricted sample of 
bright, red, centrally-concentrated galaxies, we continue 
to find a strong residual age-density relation, spanning 
~ 0.2 dex in < A5 >, while the metallicity-density rela- 
tion is diminished — although a significant residual re- 
lationship persists between metallicity and environment 
beyond that observed with color and luminosity. 

Another potential systematic effect to consider is the 
impact of aperture bias on our analysis. At z = 0.1, an 
SDSS fiber (3 arcsec in diameter) subtends a physical 
distance of only ~ 5.5 kpc (assuming a Hubble constant 
of Hp = 70 km s" 1 Mpc" 1 ). Thus, the IGallazzi et all 
(|2005f ) metallicities and ages are estimates for the cen- 
tral portion o f each galaxy How ever, as discussed in 
more detail by IGallazzi et alj (|2005f ), there is no apparent 
dependence of mean stellar age on galaxy concentration 
or redshift, at fixed luminosity within their catalog. The 
lack of any significant aperture bias is additionally sup- 
ported by detailed studies of early-type galaxies in the 
local Universe (e.g.. iMehlert , et al.'l l200a IWu et aHl200H: 
ISanchez-Blazquez et al.1 120071 ). which find no significant 
age gradient out to large physical radii. As a secondary 
check, we restrict our sample to those sources with an r— 
band half-light radius between 1 and 2 arcseconds as well 
as r < 17 and n r < 2.5. The residual age-density and 
metallicity-density relations for this subsample remain 
highly significant, as shown in Table O These results are 
consistent with a lack of aperture bias. 
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TABLE 2 

Statistical Strengths of Relations with A5 for Different Galaxy Subsamples 

Sample N~, lo Sio(*/yr) log 1() (Z/Z Q ) 

galaxy p slope CTslopc p s l opG cr slopc 

r < 17.5 104,745 O069 0333 O042 O039 0135 0.017 

r < 16.5 30,446 0.076 0.343 0.081 0.029 0.143 0.057 

r < 16.5. n r > 2.5 27,024 0.066 0.325 0.120 0.022 0.122 0.052 

r < 17, n r > 2.5, b/a > 0.7 29,002 0.049 0.288 0.096 0.030 0.085 0.033 

r < 17,n r > 2.5, 1" < R 50 < 2" 28,020 0.043 0.298 0.052 0.029 0.128 0.065 



Note. — We quantify the correlations between stellar age and metallicity and residual 
density (A5) for the set of galaxy samples analyzed in Figure [5] The number of galaxies 
(^galaxies) m each sample is given in the second column. In addition, we list the Spearman 
ranked correlation coefficient (p) as well as the measured slope (and its uncertainty) derived 
from linear fits to the corresponding binned data points in Figure [5] 



4.2. Stellar Mass versus Velocity Dispersion 

While the topic of assembly bias, as applied towards 
galaxy assembly versus halo assembly in this work, fo- 
cuses on the relationship between galaxy assembly time 
and clustering at fixed assembled mass (i.e., fixed stellar 
mass), it is also common and interesting to study differ- 
ences in galaxy clustering as a function of age at fixed 
halo mass or fixed total (dark plus baryonic) mass. For 
instance, man y studies of nearby samples in clusters and 
the field (e.g iThomas et al.l l2005l; iClemens et aLI [20061: 
iTraeer et aLll2008f) compare the ages of stellar popula- 
tions in early-type galaxies at fixed stellar velocity dis- 
persion (a common proxy for halo mass). 

Here, we employ the velocity dispersion measure- 
ments generated by the Princeton SDSS spectral re- 
duction pipeline (Schlegel & Buries in prep). The 
velocity dispersion, er, for each galaxy is estimated 
via the "direct-fit ting" method (jBurbidge et al.l 119611 ; 
IRix fc Whiti[l99l . in which the best-fit o~ is determined 
by comparing (in pixel space) a Gaussian-broadened stel- 
lar template to the observed galaxy spectrum and fitting 
for the reduced~x 2 minimum. 

As shown in Figure [6l velocity dispersion has a strong 
relationship with environment, similar to the observed 
luminosity- and stellar mass-environment relations (cf. 
Fig. [3]) — the more massive systems are commonly found 
in increasingly higher-density regions, while the average 
density also increases at low velocity dispersion largely 
due to lower-mass satellite galaxies, which reside in more 
massive (i.e., more strongly-clustered) halos. The ve- 
locity dispersion measurements utilized in Figure [5] are 
not aperture- corrected. However , applying standard cor- 
rections (e.g.. iDavies et all 119871 : Uorgensen et al.l fl995h 
yields no significant changes in the strengths of the ob- 
served relations with environment. 

To remove the relationship between environment and 
stellar velocity dispersion, we subtract the mean over- 
density at the velocity dispersion of each galaxy from 
the measured overdensity: 

A 5 (a) = log 10 (l + tf 5 )- < log 10 (l + <5 5 )[log 10 (a)] >, (3) 

yielding a new "residual" environment measure that 
traces the local overdensity relative to galaxies of like 
velocity dispersion, but possibly differing luminosity or 
rest-frame color. 

In Figure [51 we show the average dependence of this 
new residual environment, A5(ct), on velocity dispersion 
and stellar age for the r < 17.5 magnitude-limited sam- 



ple of red-sequence galaxies. By construction, A5 (cr) 
exhibits no dependence on velocity dispersion (cf. Ta- 
ble [3|), while a strong residual age-density relation is 
found, following the same dependence observed with 
A5; that is, we find that after removing the depen- 
dence of environment on velocity dispersion, there ex- 
ists a highly-significant relationship between luminosity- 
weighted mean stellar age and environment, such that 
galaxies with older stellar populations favor regions of 
higher density relative to galaxies of like velocity disper- 
sion. Furthermore, the strength of the age-A5(<7) rela- 
tion (as shown in Fig. [5]) exceeds that of the age-As rela- 
tion (as shown in Fig. [4]) for the same sample of red galax- 
ies (cf. Tableland Table [3]). The relationship between 
stellar age and environment at fixed velocity dispersion 
spans ~ 0.25 dex in As(cr), while the relationship be- 
tween age and environment at fixed color and luminosity 
spans only ~0.2 dex in A5. Thus, not only do we find a 
highly- significant relationship between age and environ- 
ment when controlling for the dependence of galaxy den- 
sity on velocity dispersion, but the residual age density 
relation is observed to be stronger at fixed velocity dis- 
persion than at fixed color and luminosity (i.e., at fixed 
stellar mass). 

4.3. Addressing the Age-Metallicity Degeneracy 

It has been understood for many years that the inte- 
grated colors and spectra of stellar populations depend 
on both metallicity and f ormation epoch (i.e., age) in a 
degen erate manner (e.g., I Worthevi 119941 : iFerreras et al.l 
fl999h : increasing the metallicity of a population at fixed 
age has a similar effect on integrated colors as increas- 
ing the age at fixed metallicity. The use of more age- 
specific and metallicity-specific spectral features, such 
as Balmer, Mg, and Fe absorption features, has par- 
tially broken this degeneracy (e.g., ITraeer et alj 119981 ; 
iFaber et al.lll999T ). However, errors in estimating stellar 
age (i.e., errors in measuring the strength of specific fea- 
tures such as H/5 absorption) remain strongly correlated 
with errors in metallicity. For instance, if a galaxy's age 
is overestimated, then its metallicity is correspondingly 
underestimated. 

Even given the correlation between errors in metallicity 
and errors in age, it is highly unlikely that our results are 
biased by such effects. To start wit h, there is no physica l 
mechanism that would cause the iGallazzilFaLl (j2005h 
age estimates to be off in an environment-dependent 
manner, thereby creating a bogus age-density relation. 
Still, if we assume that the Balmer lines (H/3 and US) 
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TABLE 3 

Statistical Strengths of Relations with A5(<t) for Different Galaxy 

Subsamples 



Sample 



slope o- Blopc p slope cr slop 



r < 17.5 -0.002 -0.003 — 0H5B — mm — craso — uuw 

r < 16.5 0.001 -0.008 0.010 0.079 0.396 0.036 

r < 16.5, n r > 2.5 -0.001 -0.004 0.023 0.072 0.426 0.050 

r < 17, n r > 2.5, b/a > 0.7 0.005 0.008 0.009 0.076 0.483 0.081 

r < 17. n r > 2.5,1" < R 50 < 2" 0.000 -0.017 0.019 0.073 0.467 0.041 



Note. — We quantify the correlations between velocity dispersion and stellar age 
and the residual density, A5(<r), for the set of galaxy samples analyzed in Figure IS1 
We list the Spearman ranked correlation coefficient (p) as well as the measured slope 
(and its uncertainty) derived from linear fits to the binned data points tracing the 
A5 (<t) [login (<t)] and A5 (cr) [log 10 (t)] relations, using the same discrete bins employed 
in Figure [6] 
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Fig. 6. — The dependence of mean absolute and residual environ- 
ment on velocity dispersion (top) and stellar age (bottom), for the 
sample of red-sequence galaxies with r < 17.5. The black circles 
connected by the solid black line correspond to the dependence of 
mean absolute environment, log-^l + iSs), on velocity dispersion or 
stellar age, while the red diamonds connected by the dashed red line 
give the respective relationship with residual environment, A5(cr). 
The green triangles connect by the dotted line in the bottom panel 
show the average relationship between residual environment, A5, 
and age, as shown in Figure [2] Even when substituting stellar ve- 
locity dispersion for stellar mass, the strong residual age-density 
relation persists. 

used to constrain age bv lGallazzi et all (|2005l ) are over- 
or under estimated due to poor emission-line correc- 
tions, such that a galaxy thought to be old is actually 
young, then that galaxy must also be more metal-rich, 
given the relationship between age and metallicity in stel- 
lar population models. Thus, even if errors in age mea- 
surements conspire to cause the observed residual age- 
density relation (a physically unmotivated scenario), any 
correction to the age estimates would produce a strength- 
ening of the residual metallicity-density relation. 



Overall, th er e is no reason to think that the 
IGallazzi et all (l2005h measurements arc off in an 
environment-dependent manner. However, as always, 
age and metallicity estimates are only as good as the 
theoretical models (i.e., the range of star-formation his- 
tories and stellar population parameters probed by the 
models). Still, any limitation or unprobed range of pa- 
rameter space associated with the models must be cor- 
related with environment at fixed mass in order to ex- 
plain our results. Thus, our general conclusion that dif- 
ferences exist in the typical stellar populations of nearby 
galaxies in high-density and low-density environs can- 
not be explained in terms of model deficiencies, as any 
limitations of the models (e.g., the exclusion of blue 
horizontal-branch stars) must be correlated with envi- 
ronment, which directly translates into a difference in 
stellar populations with environment. 

5. DISCUSSION 

In Sj3j we presented the relationships between (both 
absolute and residual) environment and an assortment 
of galaxy properties, including r— band luminosity, rest- 
frame g — r color, stellar mass, luminosity-weighted mean 
stellar age, and stellar metallicity. While we find no 
correlation between mean residual environment and stel- 
lar mass, we do observe significant relationships between 
< A5 > and stellar age and metallicity. We discuss the 
implications of these results in the remainder of this sec- 
tion. 

5.1. Comparison to Previous Studies 
5.1.1. Age-Density Relation 

As highlighted in |JTJ studies of red-sequence galax- 
ies in local clusters and within the field population have 
yielded discrepant results on the correlation between age 
and environment at fixed mass. On one hand, a variety 
of analyses studying early-type galaxies (i.e., E/S0 mor- 
phological types) in nearby clusters and the field have 
arrived at similar conclusions as our work, where red 
galaxies in highe r-density regions are typically older at a 
given mass (e7g iTrager et al.ll2000allThomas et al.ll2005h 
Clemens etaLl 120061 IBernardi et all l2006t iSmith et all 
20081 ). ' These studies utilized a variety of data sets from 

7 However, some of these analyses have detected evi- 
dence of galaxy assembly bias at very lo w (or sometim e s un- 
specified/undetermined) significance (e. g., [Thomas et aTJ 120051 ; 
IBernardi et~aTl l2006 : Clemens et"aT|[2006Tl . 
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high-resolution spectroscopic observations of relatively 
small numbers of cluster members to composite spectra 
constructed from early SDSS data. 

In contrast, more recent analyses have y i elded contra- 
dictory conclusions. While iThomas et al.l (|2005[ ) found 
cluster ellipticals to have older stellar populations than 
their field counterparts at the same velocity dispersion 
or ste llar mass, subsequent analysis by IThomas et al.l 
(|2007f ) concluded that there is no significant environmen- 
tal dependence to the slope or zero-point of the age-mass 
relation for the "bulk" of the nearb y elliptical population 
(e.g., at ages of > 4 Gyr). Instea d. IThomas et all (120071) 
claim that the trend detected by IThomas et al.1 (|2005[ ) is 
driven entirely by a population of "rejuvenated" galaxies 
(i.e., ellipticals with recent star-formation activity and 
thus ages < 3 Gyr) , which preferentially reside in low- 
density environs. 

Our results disagree, in part, w ith the more-recent 
conclusions of IThomas et al.l (|2007f ). In agreement with 
their work, we find evidence for a variation in the frac- 
tion of "rejunvenated" galaxies with environment (that 
is, we find a residual age-density relation at young stel- 
lar ages). However, as shown in Figure 2] and Figure [SJ 
even when restricting our sample to galaxies with ages of 
log 10 (i/yr) > 9.5 (i.e., excluding the "rejuvenated" pop- 
ulation), we detect a strong relationship between age and 
re sidual environment — in conflict with the conclusions 
of IThomas et~aTI ([200l . That is, we find that the corre- 
lation between luminosity- weighted mean stellar age and 
environment at fixed mass spans the entire range of ages 
probed by our various galaxy samples, rather than being 
an effect driven by only those galaxies with young stellar 
populations. 

To emphasize and further illustrate this point, we se- 
lect multiple subsamples according to various age selec- 
tions (e.g., t > 10 9 - 5 , 10 9 - 6 , 10 9 - 7 , and 10 9 75 yr). For 
each of these samples of older galaxies, we independently 
repeat our analysis — that is, we measure and remove 
the color-luminosity-density relation in each subsample 
and then examine the residual relationships between en- 
vironment and color, luminosity, stellar mass, age, and 
metallicity. As shown in Figure [7] for a sample of red 
galaxies selected to have log 10 (t/yr) > 9.6 and r < 17.5, 
we find no significant correlation between residual en- 
vironment, A5, and stellar mass. However, we do find 
significant residual age-density and metallicity-density 
relations, similar to those previously observed for sam- 
ples selected according to apparent brightness and mor- 
phology yet spanning the entire range of ages probed (cf. 
Figure [3]) . 

The results as presented in Figure [7] are independent 
of the age limit employed to select the galaxy subsam- 
ple; for samples with ages > 10 95 , 10 9 7 , and 10 9 75 yr, 
we find no significant variation in the residual relation- 
ship with environment for stellar mass, age, and metal- 
licity. Moreover, we find similar results for subsamples 
of galaxies selected according to a mass- dependent age 
criteri on, analogous to that employed by IThomas et al.l 
(120071) . When re moving the set of galax ies referred to as 
"rejuvenated" by IThomas et all (|2007l) (i.e., those with 
younger stellar populations), we still find a strong age- 
density relation at fixed stellar mass within the popula- 
tion of galaxies with older stel lar populations (i.e. , the 
"bulk" population according to IThomas et al.ll2007f ). 



Finally, when defining subsamples selected on mor- 
phology and apparent brightness as well as age (e.g., 
r < 17, n > 2.5, and log 10 (i/yr) > 9.7), we find trends 
between residual environment and stellar age and metal- 
licity analogous to those presented in Figurc[51 However, 
by excluding those galaxies with younger stellar popu- 
lations in combination with morphological and appar- 
ent magnitude selections, the number of sources at low 
mctallicities (e.g., log lo (Z/Z ) < —0.5) is significantly 
diminished such that the range in metallicities probed 
is limited and the resulting residual metallicity-density 
relation is only statistically significant within the metal- 
rich regime (i.e., log lo (Z/Z ) > —0.35). Still, altogether, 
this analysis illustrates that our results are i nconsistent 
with the conclusions of IThomas et alj (|2007t ). such that 
we find a relationship between stellar age and environ- 
ment at fixed mass that spans all ages probed in our 
study (i.e., the age-density trend at fixed mass is not 
solely driven by galaxies with young stellar populations 
preferentially residing in undcrdense environs). Addi- 
tional evidence for this conclusion, based on independent 
statistical tests, is presented in th e Appendix. 

Other recent analyses have, like IThomas et all (|2007f ). 
found little evidence of galaxy assembly bias at the mas- 
sive end of the red sequence based on archaeological in- 
vestigations of local galaxies. Studying a dozen early- 
type g alaxies in the Coma cluster in detail, iTrager et al.l 
(|2008f ) also find no difference between the inferred stel- 
lar ages of red-sequence galaxies in the cluster environ- 
me nt and simila r gala xies in the field. As discussed 
by ITrager et all (|2008D . the relatively young ages of 
the red galaxies in Coma are most likely attributed 
to recent star-formation episodes. However, nearly 
all of the galaxies in their sample show mean ages 
of 5-7 Gyr, with none having stellar populations as 
old as ~ 10 Gyr. Past spectroscopic observations of 
early-type systems in Coma have generally arrived at 
a broader range of stellar ages, with a non-negligible 
population of galaxies with ages of roughly 10 Gyr 
(IMoore et al.l 120011; IPoggianti et"all 120011; iMoore et all 
1200a INelan et al.M2005USanchez-Blazquez et al.1 12006ft . 



However, the ITrager et aL I (120081) data are of somewhat 
higher resolution and/or signal-to-noise ratio, facilitat- 
ing improved emission corrections and precision in the 
resulting age estimates. 

Accepting the accuracy of the ITrager et alj ([20081 ) age 
measurements, our results are reconcilable if either the 
ITrager et al.1 (|2008f ) sample is a biased tracer of the Coma 
early-type population (e.g., due to shot noise in the se- 
lection of their small sample) or if there are no (or at least 
very few) old early-type systems in the Coma cluster. If 
the latter is true, then Coma would be unusual amongst 
massive clusters in the local Universe. Coma is not a 
cool-core or cooling -flow cluster; instead it has a cool- 
ing time > 10 Gyr (|Kaastra et al.ll2004T ). Thus, recent 
star formation at a significant level within the central 
region of the cluster seems unlikely, especially across the 
entire early-type population. However, Coma is not a 
relaxed s ystem; instead, it contains significant subs truc- 
ture (e.g., IBiviano et aLlll996tlNeumann et £il.|[2003h and 
is ob served to be undergoing a merger/accre tion event 
(e.g-. lMellier et al.lll988t Neumann et al.ll200lD . As high- 
lighted bv lGerhard et aL I d20070, a large fraction (~30%) 
of galaxies in Coma are likely associated with the accre- 
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Fig. 7. — The dependence of mean residual environment, < A5 >, on absolute magnitude, rest-frame color, stellar mass, stellar age, and 
stellar metallicity for a sample of red-sequence galaxies with ages of > 10 9 ' 6 yr and r < 17.5. The points and error bars correspond to 
the mean environments and la uncertainties in the means computed in distinct bins of each galaxy property. While we find no significant 
residual trend between environment and color, luminosity, or stellar mass, we observe a strong dependence of residual environment on 
stellar age and metallicity. 



tion of the NGC 4839 subcluster. If Coma is a currently- 
forming supercluster, then perhaps recent episodes of 
star formation, associated with accretion/merger events, 
can explain t he rel atively young ages measured for the 
iTrager et al.l (|2008D sample of early-type galaxies. 

All told, t he young ages of ear ly-type galaxies in Coma 
as found by ITrager et alj (|2008f ) persists as a significant 
outlier relative to the results of other studies of local clus- 
ters, including other contempor ary work studying C oma. 
For example, recent results from lSmith et al. (2009]), em- 
ploying spe ctroscopic obs e rvatio ns of comparable quality 
to those of ITrager et al.l (|2008| ) and spanning a larger 
sample of early-type galaxies extending to lower masses, 
find a significant population of galaxies with old stellar 
popu lations (t ~ 10 Gyr) in the core of Coma. More- 
over, I Caldwell et al.l (|2003f ) find massive early-type sys- 
tems in Virgo with similarly old luminos ity-weighted 
ages ( t > 9 Gyr). While the sample of ISmith et al " 



2008 j is biased to lower masses than that of lTrager et al 



2008f) . one would expe ct this to lead to s maller, not 



larger, ages than for the ITrager et al.l (|2008f ) sample, as 
more massive systems tend to cont ai n older stellar pop- 
ulations (e.g . ITrager et all I2000bllat INelan et all [20051; 
Rogers et all 120081: iMichielsen et ail 120081 : ISmith et all 



20091: iMatkovic et al.l l2009) 



There are several other recent studies of local field and 
cluster galaxies that have found only weak statistical ev- 
idence for a relationship between st ellar age and e n viron - 
ment at fixed mass. For example, iGallazzi et al.l (|2006l ) 



studied the environmental dependence of the scatter in 
the age-stellar mass relationship using the same cata- 
log of stellar ages and metallicities as employed in this 
work. They find a relatively weak correlation between en- 
vironment and age such that galaxies outside of the most 
dense environments ten d to have younger ste llar popu- 
lations. The analysis of IGallazzi et a l. (2006), however, 
was limited in its scope due to a rather small sample size 
(< 2000 galax ies) for which env i ronme nt estimates were 
available fromfKauffma nn et al.l (|2004l ). Our analysis ex- 
pands upon this earlier work by markedly increasing the 
sample size with environment information and by utiliz- 
ing a m ore sensitive statistica l me thodology. 

While lClemens et al.l (|2006l) and lBernardi et all (f2006h 
both conclude that there are differences between the ages 
of massive ellipticals in differing environments locally, 
with systems in high-density regions being roughly 1 Gyr 
older at fixed velocity dispersion, evidence for a signifi- 
cant age-density relation at fixed mass in their samples is 
considerably weaker than that presented herein and more 
in line with the result s of iThomas et al.l (|2007f ) . For ex- 
ample, I Clemens et al.l ()2006l ) find a significant difference 
between the ages of ellipticals in high- and low-density 
environs over only a limited range in velocity dispersion 
(or specifically for only 2 of their 7 discrete bins in veloc- 
ity dispersion), with the differences (in t hose two bins) 
signifi cant at less than 2a (see Fig ure 10 oflClemens et all 
I2006h . Similarly, in the work bv iBernardi et al.l $00W . 
the apparent difference in average age between ellipticals 
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in high-density and low-density environs depends on the 
particular Balmer absorption line (H/3 versus Hyp) em- 
ployed to constrain stellar age, with ages based on H/3 
showing no significant difference as a function of environ- 
ment and with H7F yielding age variatio ns significant at 
the a pproximately 2a level (see Table 9 of lBernardi et al.l 
[200^ . 

T he lower significance a t which IBernardi et alj (|2006f ) 
and iClemens et al.l (120061 ) find a trend between stellar 
age and environment at fixed mass likely results from 
the smaller sample sizes, alternate methods of measur- 
ing environment, and different statistical methods em- 
ployed in the analysis. For example, by selecting sam- 
ples drawn from the extrem es of their measured envi- 
ronment distri b ution s, both IBernardi et alj (|2006f ) and 
IClemens et all (f2006h are limited to final sample sizes 
of < 90 00 and < 4000 galax ies, respectively. On top 
of that, IClemens et al.l (|2006l ) further reduce the statis- 
tical power of their sample by dividing into course en- 
vironment bins and then subdividing i nto bins of ve- 
locity dispersio n (e.g., see Fig ure 10 of IClemens et alj 
l2006f) . While IClemens et alT (120071 uses an environ- 
ment measure similar to ours. IBernardi et al.l (|2006[ ) use 
the 3-dimensional distance to the 10 th -nearest neighbor 
with M r < —21.5 as an environment estimator. By 
restricting the tracer population to relatively luminous 
galaxies, while also using a 10 th - vers us 5 th -nearest 



neighb or distance (as used in this work) . IBernardi et alJ 
(|2006f ) are tracing environment on much larger scales, 
perhaps missing the dependence of galaxy properties on 
the more local galaxy density. Moreover, measuring 
distances in 3-dimensions versus in projection leads to 
greater sensitivity to the effects of redshift-space distor- 
tions, which smear out th e galaxy distributio n in envi- 
ronments such as clusters (jCooper et al.ll200"5[) . Finally, 
alon g with using the dista nce to the 10 th -nearest neigh- 
bor, [BernaxdTeTal] (|2006| ) use the distance to the near- 
est galaxy cluster as a means for identifying high-density 
environments, thereby significantly biasing their sample 
against group-like environs. Altogether, these effects 
likely lead to the marginally-significant trends between 
stellar age and envir onme nt at fixed mass pres ented by 
IBernardi et "all (|2006f ) and lClemens et all (|2006f) . In con- 
trast, our methodology, which utilizes the entire galaxy 
sample (versus only studying galaxies in extreme envi- 
rons), yields a strong age-density relation at fixed mass. 

5.1.2. Metallicity-Density Relation 

In addition to a strong residual age-density rela- 
tion, we also find a relatively weak residual metallicity- 
density relation, where more metal-rich galaxies favor 
regions of higher galaxy density relative to galaxies of 
like stellar mass. However, the significance of this 
metallicity-density relation decreases as the sample is re- 
stricted to more centrally-concentrated (i.e., more bulge- 
dominated) or less elongated galaxies (cf. HA.1\\ , such that 
a weak (though statistically significant) correlation is ob- 
served between mean residual environment, < A5 >, and 
stellar metallicity for the bulk of the galaxy samples stud- 
ied here. 

Previous analyses of early-type galaxies within lo- 
cal clusters have generally found no evidence for a 
metallicity-density relation at fixed mass, in contrast 



to our results. For i nstance, studies of l i ne indices in 
nearb y clusters (e.g., I Guzman et al. I [19921 ; iCarter et al.l 
120021) find a significant dependence of the Mg2-<7 rela- 
tion on cluster-centric distance, where early-type galax- 
ies near the core exhibit stronger Mg2 absorption than 
their counterp a rts in the outskirts of the cluster. As 
IGuzman et alj (|1992f ) only measure one spectral index, 
they arc unable to disentangle age and metallicity; al- 
though they do conclude that a dependence of stellar 
age on cluster radius might be responsible for the ob- 
served variations in Mg2 absorption strength (sec also 
IBernardi et al.lll998f). A similar conclusion is reached 
bv iTerlevich et al.l (|200lT ) studying the U - V color of 
galaxies in the Coma cluster; color, however, is similarly 
unable to distinguish variations in age from variations in 
metallicity. 

Folding in additional s pectral indic e s to separate 
metallicity and age effects, ICarter et al.l (|2002t see also 
iKuntschner fc Daviesl fl998f ) infer a significant relation- 
ship between metallicity and radius in Coma. However, 
this metallicity-density relation does not take into ac- 
count the relationship between mass (or luminosity) and 
radius, minimizing its usefulne ss in a comparison to our 
analysis. More recent work bv lSmith et alJ ([20 06). mea- 
suring a larger number of indices (12) for a significantly 
larger sample of galaxies (^3000) in 94 nearby clusters, 
finds several significant relationships between spectral in- 
dices and cluster-centric radius at a given velocity disper- 
sion. Using these indices to distinguish between age and 
metallicity effects, they conclude that there is a signifi- 
cant age gradient in local clusters, such that systems in 
cluster cores have older stellar populations, while there 
is no significant m etallicity gradient detected (see also 
ISpolaor et afll2007)h . 

Other recent analyses spanning a broad range of envi- 
ronments in the local Universe have also found little evi- 
dence for a significant relationship between stellar metal- 
licity and local galaxy density a t fixed m ass. For exam- 
ple, b oth IClemens et al.l (120061 s ee also IClemens et al.1 
12001 and IBernardi et al l (j206H . find no significant 
metallicity-density trend at fixed velocity dispersion 
for nearby elliptica l and lenticular syst ems (see also 
iGallazzi et al.ll2006f ). Furthermore, while iThomas et al] 
(12005ft conclude that massive early-type galaxies in local 
clusters are, on average, ~0. 05-0.1 dex more metal-poor 
than their counterparts in the field, the significance of 
this result is unsubstantiated statistically, with no un- 
certainties given for the coefficients of the linear fits to 
the mean Z-a relations in the hi gh- and low-densit y 
galaxy samples (see Equation 1 of lThomas et aLlfeOOSl ). 
Considering the small size of th e galaxy sample (onl y 
124 galaxies in total) employed by lThomas et al.l (|2005h . 
the uncertainties on the parameters of the linear fits is 
likely to exceed a few percent, which would make the 
small measured differences in the mean Z-a relations in 
the t wo environment regim es statistically insignificant. 
Thus, IThomas et all ([2005), like many other authors, 
find no evidence for a metallicity-density relation within 
the nearby early-type galaxy population. 

Not only do we find a significant metallicity-density re- 
lation at fixed mass for our most inclusive color-selected 
sample, but also for our morphologically-selected sub- 
samples. While this trend is weaker than the age-density 
relation at fixed mass, it is significant nonetheless. As 
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discussed in the preceding section ( <j4.3[) . a metallicity- 
density relation is unlikely to result from correlated er- 
rors associated with the degeneracy between constraints 
on stellar age and metallicity. Moreover, when compar- 
ing apples to apples, the uncertainties in measured galaxy 
properties are generally large enough that it is far easier 
to dilute an underlying environment trend, such that it 
goes undetected, than it is to falsely detect one. Thus, 
for many previous studies of local early-type galaxies, it 
is likely that the same errors or limitations that smeared 
out the age-density relation in the galaxy sample also 
precipitated a non-detection of the metallicity-density 
relation. 

In contrast to our results and to the majority of re- 
sults fro m the aforemention e d stud ies of local early-type 
galax ies, IKuntschner et al.l (|2002t see also Ros e et al.l 
119941 ) find that their sample of nearby field ellipticals 
are more metal-rich (by ~ 0.2 dex) in comparison to 
members of the For nax cluster. This discr epancy be- 
tween the findings of IKuntschner et aTl (120021) and thos e 
of our work (as well as those of iBernardi etafl (120061) . 
iClemens et all p006l ). among others) can potentially be 
understood in terms of an age metallicity degeneracy. 
As discussed in more detail in §4.3| errors in estimat- 
ing stellar age are strongly correlated with errors in 
metallicity, such that if a galaxy's age is overestimated, 
then its metallicity will correspondingly be underesti- 
ma ted. Thus, the m e tallicit y-density relation reported 
by IKuntschner et al.l (120021) could be reconciled with 
other observations, if the IKuntschner aae 
density relation is correspondingly stronger than that 
found in other analyses. The plausibility of this scenario 
is supported by the larger difference in age between mas- 
sive ellipticals in field an d cluster environments given by 
IKuntschner et all (|2002l, ~ 2-3 Gyr) versus that found 
bv IBernardi e laTl (12001 ~1 Gyr). 8 

The analysis bv IKuntschner et al.l (|2002l ) suffers even 
more from the problem of small sample size ( < 30 galax- 
ies in total). In addition, there was a strong differ- 
ence in the morphological composition of their field and 
cluster samples, with the field population containing a 
lar ger percentage of SOs (v ersus ellipticals). However, 
the IKuntschner et~a l. (2002) study compared the ages of 
galaxies in different environments at fixed £?-band abso- 
lute magnitude and not at fixed mass (i.e., not at fixed 
stellar mass or velocity dispersion). Thus, their field sam- 
ple may have been increasingly dominated by younger 
galaxies at lower masses, which appeared bright in the 
-B-band due to recent star-formation activity. Given the 
observed correlation between mass and metallicity (inde- 
pendent of environment), such a bias in the sample would 
work to decrease the strength of the metallic ity-density 
relation detected bv IKuntschner et aTl (|2002t ). 

Recognizing the correlation between stellar metallicity 
and age observed for local massive early-type systems, 
where galaxies of a given mass with older stellar popula- 
tions tend be more metal-poo r (e.g. JGallazzi et alJl2005t 
iSanchc z-Blazo uez et al.ll2006D . it might be expected that 
an age-density relation at fixed mass would naturally be 
accompanied by some sort of metallicity-density relation, 

8 Note that Table 9 of lBernardi et al.l 120061) lists an age differ- 
ence between high- and low— density environs that is more in line 
with ~0.6 Gyr. 



like that observed bv IKuntschner et all (|2002t ). However, 
there is large scatter in the age-metallicity relation, es- 
pecially at lower masses. Furthermore, stellar age and 
metallicity measure different aspects of a galaxy's for- 
mation history; while age traces the time since the bulk 
of the stars were formed (or in the case of luminosity- 
weighted mean stellar age, the time since the last sig- 
nificant episode of star formation), metallicity reflects 
the balance of feedback and accretion processes in con- 
cert with the star-formation efficiency (that is, the effi- 
ciency with which gas is turned into stars). For example, 
two galaxies forming from pristine gas reservoirs and fol- 
lowing identical formation histories, except for an offset 
in time, will contain stellar populations of comparable 
metallicity, but distinct ages. Our results suggest that 
not only did red-sequence galaxies in high-density en- 
vironments form earlier than their counterparts in low- 
density regions, but [1] feedback and/or accretion pro- 
cesses are modulated by the local environment such that 
the high-density regime favors greater metal abundances 
and/or [2] the typical star-formation efficiency of galax- 
ies in high density regions or at higher redshift exceeds 
that of galaxies in low-density environs or at lower red- 
shift. 

Some evidence for the former scenario has been ob- 
ser yed in a stud y of sta r-forming galaxies at z ~ 0.1 
by iCooper et al.l (|2008aD . Using the gas-phase o xygen 
abundance measurements of lTremonti et al. (2004), they 
show that local environment is correlated with a non- 
negligible portion of the scatter in the mass-metallicity 
relation, such that galaxies in higher-density environs 
at a given mass tend t o be more metal-rich . As dis- 
cussed in more detail bv lCooper et al . (2008a), there are 
many physical scenarios that might lead to galaxies hav- 
ing a higher metal abundance in high-density regions 
versus in low-density regions. For example, processes 
such as ram-pressure stripping or ga laxy harassment 
(jGunn fc Gottlll972tlM"oore et al.lll996fl , which only op- 
erate in extreme high-density environs, could remove the 
outer portion (i.e., the most metal-poor portion) of a 
galaxy's gaseous halo, thereby inflating the relative metal 
abundance of gas reservoirs within galaxies in rich groups 
and clusters. Similarly, outflows and winds driven from 
galaxies in high-density environments might be less ca- 
pable of expelling enriched gas versus their counterparts 
in underdense regions. 

There is also evidence to support the possible higher 
star-formation efficiencies of star-forming galaxies at 
higher redshift relative to local star-forming galax- 
ies. Fo r example, studies of sub-millimeter galaxies 
(SMGs, iChapman et al.l 120051 and references therein) 
as well as luminous and u ltra-luminous infrared galax- 
ies (LIRGs and ULIRG s, lAaronson fc Olszewskil 1l98l 
iSanders fc Mirabelll996f ) at z > 1 typically find higher 
surface densities of star formation for a given surface den- 
sity o f molecular gas dBouche et aTl l2007t iDaddi et al.l 
l2008t iTacconi et al.l I2009D . Such a correlation between 
star-formation efficiency and redshift in concert with a 
galaxy assembly bias where galaxies in high-density en- 
virons formed their stars earlier (i.e., at higher redshift) 
than their counterparts in low-density environs would 
yield a corresponding correlation between environment 
and star-formation efficiency such that galaxies in over- 
dense regions are more efficient at converting gas into 
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stars (i.e., have higher stellar metallicity) . 

5.2. The Impact of Mergers 

The stellar populations of galaxies are not entirely 
formed via in situ star formation. Rather, much of a 
galaxy's stellar mass is assembled via accretion events 
(i.e., mergers). In addition to increasing the stellar mass 
of a galaxy, mergers can also have a substantial im- 
pact on the age of the resulting system. While gas-rich 
(or "wet") mergers will tend to induce bursts of star 
formation thereby significantly decreasing the inferred 
luminosity-weighted mean stellar age of the system, dis- 
sipationless (or "dry" ) merging of two systems with older 
stellar populations will have a less severe impact on the 
inferred age of the merger product, yielding a passively- 
evolving descendant with a stellar age that is simply the 
luminosity-weighted average of the antecedents' ages. 

In addition, mergers tend to occur in regions 
of hig her galaxy de n sity, su ch as galaxy groups 
(e.g.. iCavaliere et al.l 119921; iMcIntosh et all 120081; 
Hester fc Tasitsiomil 120091 : lldeue et all 120091 iLin et al.1 
20091) . Given this correlation with environment, it 
is important to understand the potential impact of 
hierarchical assembly on our results. Since mergers 
increase the stellar mass of the system and since galaxies 
with smaller stellar masses tend to have y ounger stellar 
popu lations (e.g., iGallazzi et al.l 120051 ; I Graves et al.1 
120091 ) . galaxies that have recently undergone a merger 
are likely to be biased towards lower ages relative 
to other galaxies of like stellar mass. In the case of 
gas-rich mergers, this effect is further accentuated by 
the decrease in age associated with the merger-induced 
burst of star formation. Thus, mergers would tend to 
decrease the average age of galaxies at a given stellar 
mass in the high-density regime, thereby decreasing the 
amplitude of the observed age-density relation at fixed 
mass. In all, the impact of mergers would be to smear 
out the residual age-density relation shown in Figure |4j 

5.3. The Evolution of Post-Starburst or 
Post-Quenching Galaxies 

Our analysis of red galaxies in the local Universe shows 
a clear relationship between stellar age and environment, 
where galaxies with older stellar populations tend to re- 
side in higher-density environments relative to younger 
galaxies of like stellar mass. Such a correlation between 
stellar age and environment could naturally explain the 
evolution in the clustering of post-starburst (otherwise 
known as K +A or post-quen ching) galaxies at z < 1, as 
observed bv lYan et all (|2009h . 

Using d a ta dra wn from the SDSS and from DEEP2, 
lYan et afl (|2009f ) studie d the distribution of e nviron- 
ments for K+A galaxies (|Dressler &; Guniilll983l ) at z ~ 
0.1 and at z ~ 0.8, relative to the red and blue galaxy 
populations at those redshifts. They find that post- 
starburst galaxies at low redshift have an environment 
distribution similar to that of blue galaxies, favoring re- 
gions less dense than those typica lly inhabited by red 
galaxies (see also iHogg et al.M2006f ): however, at higher 
redshift, they find post-starburst systems favor environ- 
ments more similar to that of red galaxies. 

If post-starburst (or K+A) galaxies have their star- 
formation quenched in the same type of environment 



(e.g., a density level corresponding to a galaxy group) 
at all redshifts less than z ~ 1, then this evolution in 
the environment distribution of post-starbursts relative 
to that of the red galaxy population from z ~ 0.8 to 
z ~ 0.1 can naturally be explained in terms of galaxy 
assembly bias. On average, a red galaxy at low red- 
shift will have had a longer amount of time since having 
its star formation quenched; thus, the environment of a 
low-z red galaxy will have had more time to increase in 
density relative to that in which the quenching occurred 
(i.e., the environment in which the K+A phase of evolu- 
tion occurs). In this scenario, the time since quenching 
(or time since assembly — thus, stellar age) for a galaxy 
on the red sequence is naturally correlated with the lo- 
cal galaxy density or environment. Our results support 
this potential explanation for the evolution of the envi- 
ronment distribution of post-starburst galaxies at z < 1. 
The observed correlation between stellar age and envi- 
ronment at fixed stellar mass supports (although does 
not prove) the picture in which K+A galaxies arc found 
in the same type of environment at z < 1. 

5.4. Implications for the Environmental Dependence of 
the Type la SN Rate 

The typical age of a stellar population has been shown 
to be strongly connected to the type la supernova (SN) 
rate, with young star-forming galaxies having much 
higher SN la rates than quiescent gala xies with older 
stellar populations (jSullivan et al.ll2006| ). The increase 
in the SN la rate with star-formation rate is commonly 
attributed to the existence of two types of la events: 
"prompt" SNc la, which likely result from the evolution 
of somewhat massive stars and are thus correlated with 
star formation, and "delayed" SNe la, which are thought 
to be the evolutionary outcome of less massive stars and 
therefore connected to the underlying stellar mass of the 
galaxy. _ 

Within quiescent systems, the delayed component of 
the la rate dominates. However, even for systems of 
like mass and lacking ongoing star formation, the SN la 
rate is still expected to depend on the age of the stellar 
population; some theoretical predictions for both single- 
degenerate and double-degenerate progenitor scenarios 
suggest that the SN la rate for a single aged stellar popu- 
lation should increase with time, reaching a peak around 
an age of ~ 10 s yr, with a sha r p decline towards late r 
times (jGreggio &: Renzhiil 119831: lYungelson et al1l!994h . 
A more recent analysis bv lGreggid (|2005h slightly revises 
this general picture, finding a sharp reduction in the type 
la rate for both single-degenerate and double-degenerate 
progenitors beyond an age of ~ 1 Gyr. 

Taking as given this theoretical connection between 
stellar age and type la rate, our results regarding the 
correlation between age and environment at fixed stel- 
lar mass along the red sequence suggest that the type 
la rate in early-type galaxies should vary with environ- 
ment. In particular, our study indicates that the SN 
la rate within quiescent early-type galaxies should be 
lower in high-density environs relative to that found in 
low density environs. The bulk of our samples have typ- 
ical luminosity-weighted mean stellar ages greater than 
1 Gyr, placing these populations beyond the expected 
peak in the type la rate and into a regime where the rate 
should be in sharp decline. 
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Currently, the evidence is mixed regarding a poten- 
tial correlation between the SN la rate (per unit mass) 
and environmen t in th e local Universe. Most notably, 
iMannucci et al.l (|2008l ) concluded that the type la rate 
is more than 3 times higher in cluster ellipticals relative 
to field ellipticals, using a sample of 11 cluster and 5 
field type la events at z < 0.04, which would run counter 
to the theoretical trend. However, the uncertainty of 
these environment-dependent rates is quite large. In 
contrast, a recent analysis of a larger sample of type la 
su pernovae drawn from the SDSS-II Supernova Survey 
bv lCooper et alj (|2009f ) found no significant relationship 
between local galaxy environment and the likelihood of 
hosting a type la event within the red galaxy population. 
Overall, such studies of the environments of type la su- 
pernovae are currently li mited by small numbers. Results 
from future surveys (e.g.. lSharon et ani2007l ; [Sand et al.l 
12001 will hopefully help in answering such outstanding 
questions. 

6. SUMMARY 

Using the measurements of luminosity- weighted mean 
stellar age and stellar metallicity from iGallazzi et alJ 
(120051) and the loc al galaxy overdensity estimates of 
iCooper et al.1 (|2009T ) . we study the relationship between 
galaxy properties and environment for various samples 
of galaxies on the red sequence. In contrast to previ- 
ous studies in the local Universe, we do not solely focus 
on the early-type galaxy population or on galaxies se- 
lected to reside in particular environments (e.g., specific 
clusters such as Virgo or Coma); instead, we define vari- 
ous galaxy samples drawn from the SDSS, controlling for 
correlations between environment and color, morphology, 
mass, etc., while probing a broad and continuous range 
of environments. Our principal results are as follows: 

1. After removing the mean color-luminosity- 
environment relation on the red sequence, we find 
a strong residual relationship between environment 
and stellar age (cf. Figure [4]), such that galaxies 
with older stellar populations favor regions of 
higher galaxy overdensity relative to galaxies of 
like color and luminosity (i.e., like stellar mass). 

2. When using subsamples restricted according to 
morphological, brightness, and color criteria, we 
find this residual age density relation persists with 
no significant change in strength (cf. Figure O and 
Tabled]). 

3. Using stellar velocity dispersion to trace mass (in 
lieu of stellar mass), we find that the age-density 
relation at fixed mass persists (cf. Figure [6]) and is 
even more statistically significant (cf. Table Q] and 
TableEJ). 

4. We find a significant residual correlation between 
environment and stellar metallicity at fixed color 
and luminosity (cf. Figure [4]). While a strong resid- 
ual Z-density relation is observed for our most in- 
clusive sample of red-sequence galaxies, enforcing 
selection criteria based on concentration, axis ra- 
tio, and/or brightness yields samples that exhibit 
weaker (though still statistically significant) resid- 
ual metallicity-density relations (cf. Figure [5] and 
Table H. 



5. Our results are in confli ct with the gen e ral co nclu- 
sions of rece n t wor k by iThomas et alj (|2007t ) and 
iTrager et al.1 (|2008t ). which both indicate a lack of 
correlation between age a nd environment at fixed 
mass. While the results of ITrager etHI ([2001 may 
be due to the peculiarities of the Coma cluster (rel- 
ative to other l ocal clusters) a n d/or sample selec- 
tion effects, the I Thomas et al.l (120071 ) findings are 
clearly not supported by our analyses. While we 
find evidence for a variation in the fraction of "re- 
juvenated" (or young) galaxies wit h environment 
(in a greement with the results of IThomas et al.l 
l2007h . we also find a correlation between age and 
environment at fixed mass even for the "bulk" (or 
old) earl y-type galaxy p o pulat ion — counter to the 
claims of IThomas et al.l (120071 ) . 

6. In agreement with several previous studies of early- 
type systems in the local Universe and with a mul- 
titude of studies at intermediate redshift, we find 
significant evidence for assembly bias on the red 
sequence, such that red, early-type galaxies that 
formed earlier in the history of the Universe are 
more strongly clustered today than their counter- 
parts (of equal mass) that formed later. 

As discussed in EJXJ evidence for an assembly bias in the 
galaxy formation process has commonly been explored in 
two distinct manners: an evolutionary approach, which 
studies galaxies over a range of redshift and directly in- 
fers evolution in the galaxy properties, and an archae- 
ological or paleontological approach, which studies the 
stellar populations of nearby galaxies in detail with the 
aim of deriving the evolutionary history from the cur- 
rent fossil record. As evidenced by the work presented 
herein, the results from these two methods are now arriv- 
ing at a relatively consistent picture. Studies of galaxies 
at z < 1 have shown that the population of red (or early- 
type) galaxies was preferentially built-up in overdense 
environments, such that red-sequence galaxies in dense 
environments have generally assembled their stellar mass 



in less dense regions (e.g 


..iTreu et al.ll2005t 


Bundv et al. 


2006 


: ICooper et al.ll2007 


:IGerke et al.H2007l: 


Capak et al. 


2007 


). In addition, studies of galaxy populations at yet 



higher redshifts (z ~ 2-3) have detected populations 
of massive gal axies with passively-evolving stellar pop- 



ulations 



iPapovich et all l2006t iStutz et all I2008D t hat are very 
strongly clu stered (e.g.. iDaddi et al. [ i2lm lOuadrT et al.l 
[20071 120081) . 

Our results, which apply the aforementioned archae- 
ological method, show a general agreement with these 
results from studies at higher redshifts, supporting a pic- 
ture of assembly bias in the local galaxy population. We 
find that at a given mass on the red sequence galaxies 
with older stellar populations are found in more over- 
dense environments (i.e., are more strongly clustered) 
than their counterparts with younger stellar populations. 
In contrast to some previous studies, we find that this 
correlation between age and environment at fixed mass 
is not solely driven by a minority population of galaxies 
that have experienced recent episodes of star formation. 
Instead, we find a strong residual age-density relation 



Labbe et all [20051: ILonghetti et all [20051: 
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even among those red-sequence galaxies dominated by 
older stellar populations (i.e., with ages >5 Gyr). 

In the future, we look to quantify the variation in the 
age of stellar populations as a function of environment, 
with the goal of constraining theoretical models of galaxy 
formation and evolution. As shown in the Appendix, pre- 
liminary analysis of only those galaxies with older stel- 
lar populations (i.e., log 10 (f/yr) > 9.75) clearly demon- 
strates that there is a difference of at least 0.25 Gyr 
between those systems residing in high-density environs 
and those in low-density reg ions. Some recent analyses 
at intermediate redshift ( e.g.. Ivan Dokkum fe Ellisll2003l: 
van der Wei et all [20051 : | van Dokkum fc van der Marell 
200l " suggest that while early-type galaxies in clusters 
are older than their field counterparts, the difference in 
stellar age is less than would be expected from current 
generations of galaxy formation models. More quanti- 
tative analysis at low and intermediate redshift will be 
required to definitively measure the relationship between 
stellar population parameters and environment at z < 1 
and thereby further constrain the theoretical models. 
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APPENDIX 

In this Appendix, we investigate further the impact of galaxies with relatively young stell ar populations on the 
measured relationship between age and environment at fixed stellar mass. As discussed in g5J IThomas et "all (|2007f ) 
conclude that the observed dependence of mean stellar age on environment at fixed stellar mass or velocity dispersion 
within the local early-type population is entirely driven by "rejuvenated" galaxies, that i s, ear ly-type galaxies with 
relatively young stellar populations (log 10 (f/yr) < 9.6 from Figure 1 of IThomas et al.l I2007T) . For the "bulk" of 
the galaxy population (i.e., the quiescent early-types with older stellar populations) . IThomas et "all (|2007f ) claim no 
correlation between age and environment at fixed mass. 

To determine if there is an age-density relation at fixed mass within the "bulk" population, we define a subsample 
of red-sequence galaxies at r < 17, with early-type morphologies [n r > 2.5), and with stellar ages above the limit 
log 10 (t/yr) > 9.6, so as to exclude the "rejuvenated" population. From this parent population, we then select those 
galaxies with stellar masses of 9.8 < log 10 (M»//i _2 M Q ) < 11 (using the KCORRECT mass estimates) and residing in 
high-density environs such that log 10 (l + 65) > 1.2. This sample of 3291 galaxies represents the extreme high-density 
tail of the environment distribution (comprising < 8% of the r < 17.5 red sequence population within this mass and 
age regime). As a comparison sample, we also select those galaxies within the same apparent magnitude, morphology, 
stellar mass and age limits and with measured local overdensities below the cosmic mean (i.e., log 10 (l + 65) < 0), or 
the least overdense ~40% of the sample. Now, from this relatively large set of galaxies in low-density environments, 
we randomly draw galaxies so as to match the stellar masses and redshifts of the high-density subsample. That is, for 
each galaxy in the high-density subsample, we randomly draw a member of the low-density subsample from within 
a 2-dimensional window of A log 10 (M*//i _2 MQ) 2 + Az 2 < 0.002, centered on the stellar mass and redshift of the 
particular galaxy in the high-density subsample. Some galaxies are duplicated in this low-density comparison sample; 
however, the relatively large size of the parent low-density subsample ensures that duplication is minimal such that 
> 90% of the comparison sample is unique, with no individual galaxy included more than 3 times. 

For these two samples, matched in stellar mass and redshift but drawn from distinct portions of the environment 
distribution, we then compare the distributions of stellar masses, ages, and metallicities. As shown in Figure [51 the 
stellar mass distributions for the high-density and low-density samples are in excellent agreement, as expected given 
the construction of the populations. However, the distributions of stellar ages for the two samples are quite distinct, 
with the low-density population skewed to younger ages relative to their high-density counterparts. Similarly, the 
galaxies in low-density regions tend to be more metal-poor. 

Applying non-parametric statistical tests to the data support the perceived differences in the binned age and 
meta llicity distributions show n in Figure [SJ Specifically, we use the one-sided Wil coxon-Mann-Whitney (WMW) U 
test (IMann fc Whitnevi 11947ft and the two-sided Kolmogorov-Smirnov (KS) test (press et all 119921 ; IWall fc Jenkins! 
12003ft in our analysis. The result of each test is a P-value: the probability that a value of the WMW U or KS statistic 
equal to the observed value or more extreme would be obtained, if some "null" hypothesis holds. Results with a 
P-value below 0.05 (roughly corresponding to 2er for a Gaussian distribution) are considered to be significant, while 
w e term P-values less than 0.01 as highly significant. For more details regarding the WMW U and KS statistics, refer 
to lCooper efaLl (|2009ft . 

Performing a one-sided WMW U test on the age distributions for our high-density and low-density subsamples, we 
find that the galaxies in low-density regions do have younger ages, with a P-value <C 0.01. Thus, there is a much 
less than 1% chance that we would observe a difference this strong if both samples were drawn from the same parent 
distribution. For the KS test, we arrive at a similar result, also finding a P-value <C 0.01 when comparing the age 
distributions for the galaxies in overdense and underdense environs. Thus, when limiting the sample under study 
to galaxies with older stellar populations (i.e., the "bulk" population), we find that red-sequence galaxies residing 
in low-density environs are younger than galaxies of like stellar mass in high-density regions. Applying these same 
tests to the metallicity distributions, we again find highly significant results, with P-values <C 0.01, confirming that 
galaxies in low-density regions are metal-poor relative to their counterparts with the same stellar mass in high-density 
environments. 

Repeating the analysis above, while increasing or decreasing the stringency of the morphological-, r-band 
magnitude-, and age-limits applied to the parent galaxy sample, we find that the WMW U and KS tests yield 
similar results (cf. Tabled]). F° r example, when limiting the parent population to those galaxies with increasingly 
older stellar populations (i.e., r < 17, n r > 2.5, and t > 10 9 75 yr), the distributions of ages for the high-density 
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TABLE 4 
Results of Statistical Tests 



Sample 


■^overdense 


lo glO 


(M.) 
Pks 


lo gl0 
AvMW 


(t/yr) 

Pks 


log lo (Z/Z ) 

Avmw ^ks 


r < 17.5, t > 10 9B 
r < 17.5, t > 10 9 ' 75 


6016 
5266 


0.40 
0.37 


0.87 
0.84 


0.00 
0.00 


0.000 
0.000 


0.000 0.000 
0.000 0.000 


r < I7,n r > 2.5, t > 10 M t> 
r < 17, n r > 2.5, t > 10 9 75 


3291 
3016 


0.47 
0.39 


0.95 
0.92 


0.00 
0.00 


0.000 
0.000 


0.000 0.000 
0.000 0.000 


r < 17, n r > 2.5, b/a > 0.7, t > 10 U B 
r < 17, n r > 2.5, b/a > 0.7, t > 10 9 ' 75 


2016 
1901 


0.40 
0.46 


0.92 
0.96 


0.00 
0.00 


0.00 
0.00 


0.000 0.000 
0.000 0.000 





Note. — We tabulate the P-values, Pwmw an d ^KS; from comparing the stellar mass, age, and 
mctallicity values for the low-density and high— density subsamples drawn from the listed parent samples, 
using the Wilcoxon-Mann— Whitney (WMW) U test and the two-sided Kolmogorov-Smirnov (KS) test. 
As discussed in the text, smaller values indicate a lower probability that the observed differences in the 
samples will occur by chance if they are selected from the same underlying parent distribution. Note 
that, as a one-sided test, the P— values for the WMW U test have a a maximum value of 0.5. For each 
sample, the number of galaxies in the overdense subsample (N oveI< ^ ense ) is listed as well. 

and low-density samples are distinct at a highly-significant level (again, P <C 0.01). As shown in Figure [51 these 
log 10 (f/yr) > 9.75 samples have age distributions well fit by Gaussians, lacking any tail towards lower stellar ages and 
thus sample the "bulk" population exclusively. 

The results i llustr ated and enumerated in Figure [5] and Table SJ which are inconsistent with the conclusions of 
iThomas et al.l (|2007| ) , illustrate that there is a significant bias in the assembly of local early-type galaxies with older 
stellar populations such that galaxies in low-density regions formed later than their counterparts in high-density 
environments. Computing the difference between the mean and median of the age distributions for the low- and 
high-density subsamples with t > 10 9 75 yr, we find an offset of roughly 0.25 Gyr. This, however, is a lower limit 
to the age difference between early-type galaxies in underdense and overdense environments. Measurements of local 
environment arc inherently noisy (relative to other measures of galaxy properties like rest-frame colors, luminosities, 
and even stellar masses), causing correlations with the local galaxy density to be smeared out. Also, by using a broad 
range in log 10 (l + S 5 ) over which to select our low-density subsample, we reduced the dynamic range of any intrinsic 
age-environment correlation. 
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Fig. 8. — The distributions of stellar masses (left), ages (center), and metallicities (right) for red-sequence galaxies in high-density 
environments (solid black line) and for the comparison sample of galaxies in low— density regions, selected to have matching stellar masses 
and redshifts (red dashed line). The two environment-specific subsamples are selected from a parent population with 9.8 < log 10 (M«) < 11, 
r < 17, n r > 2.5, and log 10 (i/yr) > 9.6 (top) or log 10 (t/yr) > 9.75 (bottom). Even when excluding those galaxies with young stellar 
populations, there remains strong evidence for assembly bias on the red sequence. 



